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Context. MHD turbulence with Hall effect. 

Aims. Study how Hall effect modifies the quenching process of the electromotive force (e.m.f.) in Mean Field Dynamo (MFD) theories. 
Methods. We write down the evolution equations for the e.m.f. and for the large and small scale magnetic helicity, treat Hall effect as a 
perturbation and integrate the resulting equations, assuming boundary conditions such that the total divergencies vanish. 
Results. For force-free large scale magnetic fields. Hall effect acts by coupling the small scale velocity and magnetic fields. For the range of 
parameters considered, the overall effect is a stronger quenching of the e.m.f. than in standard MHD and a damping of the inverse cascade of 
magnetic helicity. 

Conclusions. In astrophysical environments characterized by the parameters considered here. Hall effect would produce an earlier quenching 
of the e.m.f. and consequently a weaker large scale magnetic field. 

Key words. Magnetohydrodynamics and plasmas - Plasma turbulence 



1. Introduction 

The origin and evolution of magnetic fields observed in all 
objects of the universe is one of the main problems in astro- 
physics. The basic physical process assumed to create them is 
a dynamo, which needs two basic ingredients, a seed field and 
an amplifying mechanism, each of them constituting at present 
an independent line of research (e.g., Grasso & Rubinstein 
nMT\ Widrow lMHI GiovanninignOS- An amplifying mecha- 
nism usually considered is the so called turbulent Mean Field 
Dynamo (MFD), as turbulence is normally present in astro- 
physical environments. In this mechanism it is assumed that 
turbulence is excited at a small scale £s and that as a conse- 
quence a magnetic field is induced at a larger scale i^. This 
theory has been a useful framework for modeling local origin 
of large scale magnetic fields in stars and galaxies. 

In the Universe there are very different astrophysical en- 
vironments: compact stars, low density and low temperature 
plasmas, accretion disks around stars and in AGN's, etc. The 
plasma in each of those ambients has a different composition 
and therefore different physical processes may be relevant: in 
low ionized plasmas as the interstellar medium, ambipolar di- 
fussion is important (Zweibel 1988 1; in the high-temperature 
intracluster gas ohmic dissipation plays a major role, and Hall 
effect can be relevant in accretion disks (Sano & Stone 2002i 
WardIe|T999 Balbus & Terquem 2001 1 as well as in the early 
universe (Tajima et al. ll992t . Turbulent dynamo operation may 



therefore be affected by the composition of the plasma: if we 
consider a plasma formed by, e.g. protons, electrons and neu- 
trals, then the different interactions among these constituents 



can be expressed as a generalized Ohm's law (Spitzer [19621 
Priest 1982). 

One of the main steps in the development of a turbulent 
mean field (or large scale) dynamo theory was the recognition 
of the pivotal role played by magnetic helicity (e.g., Pouquet et 
aI. [T976l Blackman & FieId l2001all200Tbl Brandenburg lMHt . 
In the absence of resistive dissipation, and for boundary condi- 
tions such that total divergencies vanish, this quantity is glob- 
ally conserved, independently of any assumption about the tur- 
bulent state of the system. Its evolution does not explicitly de- 
pend on the non-linear backreaction due to Lorentz force, it 
merely depends on the induction equation, providing therefore 
a strong constraint on the nonlinear evolution of the large scale 
magnetic field. 

As stated above, mean field dynamo amounts to split the 
fields into large scale mean fields Uq, Bq, Aq and small scale 
turbulent fields u, b, a. This small scale fields represent the 
"waste product" of turbulence, and they can be very intense 
in spite of their small coherence length In this theory the 
evolution equation for Bo can be cast as dBo/dt = V x 
(Uo X B() -H e - rjJo), where Jo is the mean electric current, rj 
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' In order to generate large scale fields, helical turbulence is needed. 
The generation of these small scale fields is not to be considered as a 
result of a small scale dynamo. Those dynamos require turbulence to 
be non-helical (see e.g., Zel'dovich et al .l983j . 
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the resisitivity and e = (u x b)() the turbulent electromotive 
force (e.m.f.). In the two scale approach it is assumed that e can 
be expanded in powers of the gradients of Bq in the rather gen- 
eral form Si = aij [g, Cl, Bq, ...) Boj+jjijk (g, Cl, Bq, ...) dBoj/dxt, 
where the functions Oij and rjij^ are called turbulent transport 
coefficients. They depend on the stratification g, angular ve- 
locity n and mean magnetic field Bo.They may also depend on 
correlators involving the small scale magnetic field in the form, 
for example, of small scale current helicity. 

The simplest way of calculating the turbulent transport co- 
efficients consists of linearizing the equations for the small 
scale quantities, ignoring quadratic terms that would lead to 
triple correlations in the expressions for the quadratic terms. 
In other words, the backreaction of mean field Bq on the cor- 
relation tensor of the turbulence is taken into account, while 
neglecting the effect of the small scale fields. In this way 
the electromotive force can be written as (Krause & Radler 
nM}> e = aBo - j8Jo with a ^ -(1/3)tco„-<u ■ V x u>o and 
/3 ^ {l/3)Tcorr{u^}o, with Tcorr the coiTclation time. These 
modifications of the turbulent transport coefficients have been 
calculated about thirty years ago, and the approximation is 
known as First order smoothing approximation, or FOS A (see 
also Moffat \WTn Rudiger \W7^ Parker [1955 Moffat W7^ 
Zel'dovich, Ruzmaikin & Sokoloff fT983l . 

There remains to incorporate the modifications to s that in- 
volve the small scale, fluctuating fields. They arise when cal- 
culating (u X b)o from terms involving the nonlinear terms 
and the Lorentz force in the evolution equations for b and 
u respectively. As a consequence, the a term written above 
gets renormalized in the nonlinear regime by the addition of 
a term proportional to the current helicity (b ■ (V x b) of the 
fluctuating field, which in turn is related to the magnetic he- 
licity of the small scale magnetic field. The /3 term on the 
other side is not affected by the backreaction of the small 
scale fields (Pouquet et all 1976 Subramanian & Brandenburg 
120041 Brandemburg & Subramanian 2005 ). In this way we have 
a - -(l/3)w((u ■ (V x u)>o - <b ■ (V X b)>o). 

In this paper we investigate how Hall effect modifies the 
process of quenching of ffie e.m.f. s, described in the previous 
paragraphs. For this purpose we use a closure scheme recently 
introduced by Blackman & Field (2i)02 2004) . that permits to 
study dinamically the backreaction of both the large and small 
scale fields on s. This closure, also named the "minimal r ap- 
proximation", consists in finding the evolution equation for the 
electromotive force instead of finding s itself. In it, three point 
correlations of the generated small scale fields, T, are not ne- 
glected, but their sum is assumed to be a negative multiple of 
the second order coiTelator, i.e. T - -s/t. This assumption 
produces results that are in very good agreement with numeri- 
cal simulations (Brandenburg & Subramanian 2005J. 

Hall effect is taken into account by considering a corre- 
sponding term in Ohm's law (see Spitzer [79621 Priest 19821. 
The parameter that meassures its intensity is the Hall length, 
that in alfvenic units is defined as iff - (47rp)'''^ /n^e, with p 
the mass density and the electronic numerical density. Of in- 
terest is the ratio of this length, to the Ohmic dissipation length, 

and to the scale of the flow, 4. For in < I,,, the Hall-MHD 



equations reduce to those of standard MHD, as ohmic dissipa- 
tion erases any other interaction. In several astrophysical prob- 
lems, such as accretion disks, protoplanetary disks, the early 
universe plasma and the magnetopause (Birn et al 2001 , Balbus 
& Terquem 2001 Sano & Stone E002l Tajima et al[T^2j, the 
Hall scale is larger than the Ohmic scale, but it can be smaller 
or larger than 4. 

Under the hypothesis that large scale fields are force-free 
(which is a reasonable assumption in many astrophysical en- 
vironments) we obtain evolution equations for the mean mag- 
netic field and for the large and small scale magnetic helici- 
ties, that are formally identical to those obtained in absence of 
Hall effect, provided that we redefine the electromotive force 
by bh - {Ue X b) with - u - (h{(^ x b) x b)o. This means 
that the component of eh along the mean field Bq governs the 
evolution of magnetic helicity, and in turn magnetic helicity in- 
fluences the growth of Bq ( Parker ll955l JiLl999l. We study a 
system of coupled evolution equations for eh and for the large 
and small scale magnetic helicities. As mentioned above, the 
evolution equations for the helicities are formally identical to 
the ones in absence of Hall effect. In contrast, the equation for 
Eh presents substantial differences in comparison to the stan- 
dard MHD case: (i) In the a term, proportional to Bo, the fluid 
helicity is replaced by the electronic fluid helicity and there 
appears an extra term, explicitly dependent on In that couples 
b with VLe- (ii) In the turbulent diffusion term, proportional to 
V X B(), the /3 term, the fluid kinetic energy term (m^)) is re- 
placed by (u ■ Uf.)() and there also appears a correction explicitly 
dependent on £h, that couples b with u. All these modifications 
render this term not positive definite, a fact that could produce 
a transfer of energy from small scales toward large scales or, in 
other words an inverse cascade of energy (Mininni et al. l20051 . 
(Hi) There appears a new term, proportional to V^Bo, that again 
couples the mentioned velocities. The coupling of Ue with b 
indicates that Hall effect acts by transferring energy between 
these two fields in a non-trivial way. 

In order to illustrate how Hall effect affects the quench- 
ing process of the mean field dynamo, we applied the obtained 
equations to a specific physical situation in which we consid- 
ered that turbulence of maximal kinetic helicity is excited at a 
certain scale ^ and that the large scale magnetic field is gener- 
ated at a scale ti = 5^'.,. As for the Hall effect, we considered 
{s < tn < (l and treated it as a perturbation. We find that 
for this situation the overall effect is a quenching of the e.m.f. 
stronger than in standard MHD, acompanied by a supression 
of magnetic helicity inverse cascade. Our results are in quali- 
tative agreement with recent numerical simulations performed 
by Mininni et al (I2003b> . 

As the aim of this paper is to understand conceptually how 
Hall effect acts on a MFD, we did not apply our results to a con- 
crete astrophysical object. We leave this issue for future work, 
after a deeper understanding of the mechanism is attained. The 
paper is organized as follows: in section §2 we present the 
main equations and deduce the evolution equations for large 
and small scale fields. In section §3 we deduce the dynamo 
equations, i.e., the ones for the stochastic electromotive force 
and for the large and small scale magnetic helicities. In sec- 
tion §4 we implement a two scale approximation and numeri- 
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cally integrate the system of equations and discuss the results. 
Finally in section §5 we sumarize our conclusions. 

2. Main Equations 

In Magnetohydrodynamics, Hall effect can be taken into ac- 
count through the generaUzed Ohm's law as (e.g., Spitzer .1962 
PriestEnH: 



E + U X B = — J X B + T/V X B , 



(1) 



with J = V X B, «(, the electron number density, e the modu- 
lus of the fundamental electric charge and 77 the Ohmic diffu- 
sion coefficient. We need the magnetic field induction and the 
Navier Stokes equations. We use units in which the magnetic 
field has dimensions of velocity. To simplify the calculations 
and the comparison with previous works we shall consider an 
incompressible fluid, i.e., V ■ U = 0^. This condition is ful- 
filled in several astrophysical environments. The Navier-Stokes 
equation is then written as: 



— = -P [(U ■ V)U - (B ■ V)B] - vV X (V X U) , 

at 



where P s I - VV-/V^ is the projector operator onto the sub- 
space of solutions of this equation that satisfy the condition of 
incompressibility (McComb 2003J and v the kinematic viscos- 
ity. The induction equation reads: 



5B 

— = V X {U X B 

dt 



(V X B) X B - 77V X B} 



(3) 



where we defined the Hall length. In, as iu - {^np) ' jrice. We 
also need the equation for the vector potential A. If we choose 
to work with the Coulomb gauge, i.e., V ■ A = 0, it reads: 



5A 
'dt 



P [U X B - (V X B) X B] - 77V X B , 



(4) 



where P is the previously defined projector, but now projecting 
onto the space of functions that satisfy the chosen gauge. When 
T] -0, Equation Q represents the freezing of the magnetic field 
to the electron flux. To see this, let us write U^, = U - in^ x B, 
which when substituted in eqs Q and Q transforms them in 
equations formally identical to the ones without Hall effect. 

2.1. Large and Small Scale Fields 

As we are interested in studying mean field dynamo, we split 
the fields U, B and A as U = u, B = Bq -H b and A = Aq -H a. 
Upper case and subindex denote large scale fields, i.e. vec- 
tor quantities whose value may vary in space but whose di- 
rection and sense are almost uniform or vary very smoothly. 
Technically speaking, they represent local spatial averages. 



^ This condition implies that the time that a sound signal takes to 
travel through a given distance / must be small compared to the time t 
during which the flow changes appreciably, i.e., t » l/c„ so that the 
propagation of interactions in the fluid may be regarded as instanta- 
neous (e.g., Landau & Lifshitz .l997j 



Lowercase denotes small scale, stochastic fields, i.e. fields 
whose amplitude may be large, but that have a very small co- 
herence length. We assume that any average of stochastic quan- 
tities is zero. Observe that we assumed Uq - 0, i.e. no large 
scale flows. 

2.1 .1 . Evolution equation for tine mean fields 

To derive the evolution equations for the large scale fields, we 
replace the previous decomposition into eqs. Q and @ and 
take local spatial averages that we denote as (...)o If besides 
we demand large scale fields to be force-free, we obtain: 



^ = V X + T/V^Bo , 

ot 



Mo 
dt 



—^^PsH + 77V' Ao , 



(5) 



(6) 



where with the aid of Reynolds rules (McComb l2003> we have 
interchanged derivatives with averages. We have also defined a 
Hall turbulent electromotive force as: 



(2) fiH = (u X b)o - ^//((V X b) X b)o = <u, x b)o . 



(7) 



2.1 .2. Evolution equations for the small scale fields 

The evolution equations for the small scale fields are obtained 
by replacing the decomposition of fields into global averages 
and stochastic component, into eqs. (|2}, (|3} and (0} and sub- 
stracting from them the equations for the mean fields. Thus we 
have: 



dn 

Tt 



db 

dt 

and: 



P[(b • V)Bo -H (Bo ■ V)b] -H P[(b ■ V)b - (u ■ V)u] 
vV X (V X u) , (8) 

V X (Uf X Bo) - X [(V X Bo) x b] h- V x (u^ x b) 
<V X (u, X b)>o - 77V X (V X b) , (9) 



P(u, X Bo) - {h^ [(V X Bo) X b] + P(u, X b) 

<P(u. xb))o-77Vx(Vxa) . (10) 



3. Dynamo equations 

Now that we have the complete set of evolution equations for 
large and small scale quantities, we can proceed to derive the 
evolution equations for the electromotive force and the mag- 
netic helicity. 

^ In the absence of large scale flows, i.e., if Uo = 0, we obtain from 
eq. O the following constraint: <P [(U ■ V) U]>o - <P [(b • V) b])o = 0, 
that must be satisfied in order to guarantee the vanishing of Uo for all 
time. 
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3.1. Magnetic helicity evolution equation 



+ v<V2uxb>o + T, 



(13) 



Magnetic helicity is defined as the global average, or average 
over the entire volume of A ■ B, that we denote by Hj - (A ■ 
B),.o/ (Biskamp '1997V These quantities do not vary in space, 
they depend only on time. Call H'^ = (Ao ■ Bo)ro/ and = 
{a ■ b}voi- By taking the time derivative of these quantities with 
respect to time and using eqs. 0, © and ilOi we obtain 



and: 



M 



2{eh ■ Bq),oI - 277(Bo ■ (V x Bo))™/ 



+ (V ■ [sh X Ao - 77 (V X Bo) x Ao]) 



ml : 



(11) 



dh 

— = -2(eH ■ Bo).,„, - 277«(V X b) ■ b)o).,„, 
ot 

+ <V ■ {[u, X Bo - (V X Bo) X b] X a}),,,, 
+ <V ■ {[u, X b - 77V X b] X a})yoi , 



(12) 



To deal with the operator P we followed the procedure deviced 
by Gruzinov & Diamond ( 1995) , that consists in transforming 
Fourier the equations before taking averages, and make a devel- 
opment to first order in k^/ks with ki the scale of Bo and ks the 
scale of the small scale fields. We also made some simple alge- 
braic manipulation to put the total divergencies in evidence. If 
we add up eqs il 1> and J12t we see that total magnetic helic- 
ity is conserved, except for the divergencies and the dissipative 
terms. This means that the term {sh-Ro}voI transforms magnetic 
helicity between mean and fluctuating fields. In what follows 
we consider boundary conditions such that the divergencies in 
eqs. il 1> and il2i vanish. This selection is debatable, however, 
in view of the fact that such conditions may not be quite gen- 
eral, or easily attainable in practice. Nevertheless they have two 
advantages: First, the resulting magnetic helicity is gauge in- 
variant and second, they are widely used in numerical simula- 
tions, a fact that will facilitate comparisons with those works. 
The effect of boundary conditions on the evolution and gauge 
invariance of magnetic helicity is discussed in Berger & Field 
( 19841, Ji ( 1999), Vishniac & Cho (|200L|, and Subramanian & 



Brandenburg JMHb . 

3.2. Evolution equation for 

According to its definition, eq. 0, sh is the combination of two 
terms. So we need to find evolution equations for each term and 
then join them into one equation. The derivation is sketched 
in Appendix A, eq. jA.4> . and here we quote the final result, 
namely: 



with T representing the small scale field, three-point correla- 
tions and given in Appendix A by eq. iA.3\ . The term pro- 
portional to Bo is similar to the "a" term that appears in the 
kinematic dynamo, except that now it has the electronic ki- 
netic helicity (P' term inside braces) instead of the fluid kinetic 
helicity of the ordinary dynamo. Besides this term, there is a 
current helicity term (2'"^ term inside braces, also present in 
standard MHD dynamo equation) that is due to the small scale 
magnetic field and that can be cast in terms of the small scale 
magnetic helicity. Finally there is a new term that is an explicit 
Hall modification (3'''^ term in braces), that couples the small 
scale magnetic field to the electronic velocity field. 

The term proportional to (V x Bo), named term" is also 
strongly modified: the first term turns out to be the scalar prod- 
uct of the fluid and kinetic velocities, and there appears a sec- 
ond term, explicitly dependent on £h that couples u to b. All 
these modifications made this term not positive definite any- 
more. A negative value of this coefficient represents non-local 
transfer from small scale turbulent fields to the large scale mag- 
netic field (Mininni et al I2005> . Finally there appears a new 
term, proportional to V^Bq. 

From eqs. il It and (I12t we see that the important quan- 
tity in the mean field dynamo operation is the component 
of Sh parallel to Bo (Parker 1955 1, which can be written as 
s'lf - Sh ■ Bo/IBol and whose evolution equation is then given 
by 56^5? - (5eH/5f) (Bo/IBol) + sh [5(Bo/|Bo|) Idt]. To nu- 
merically integrate the resulting equation it is more convenient 
to write Ug back in terms of u and b. The physical reason is that 
u is the velocity that can be externally excited or prescribed. 
Thus we shall work with (see ea. lA.2l of Appendix.) 



c?e" 1 

= - {<(V X b) ■ b)o - <u ■ (V X u))o - ^//<u ■ V2b)o 

+ tnib ■ W^u)o + fnii^ x b) ■ (V x u))o 

- 4(b ■ [V X v2b])o + 4((V X b) ■ v2b)o } |Bo| 

- T <v'^)o + ^h(V ■ (u X b))o 1^- 

J |i>ol 

(V^Bo) ■ Bo 



+ -{h [<u ■ b)o - ^H<(V X b) ■ b)o] 



IBol 



(14) 



where the last term represents the dissipative terms and the 
more important the three-point correlations of the generated 
small scale fields. 



dsH 
dt 



4. Solving the System 

= T [-{ue ■ (V X Ue))o + ((V X b) ■ b)o + inib ■ V^u,)o}Bo 41 Further approximations and numerical integration 



3 

- i [<u ■ u,)o + ^H<b ■ (V X u))o] (V X Bo) 



b)oV2Bo 

+ 7? [<u, X W^bh - ^h((v X V^b) X b)ol 



To numerically integrate the equations, we assume that full he- 
Ucal turbulence is excited at a certain scale ^5 smaller than the 
system's size, and that the large scale magnetic field is induced 
at a larger scale {l = 5^^, that can be the system's size. This 
assumption enables us to consider that spectra of small scale 
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quantities peak at wavenumber - 2nl{s while large scale 
quantities do so at = 2nl(i. This assumption in based on the 
work of Pouquet, Frish & Leorat ('1976), who several years ago 
showed that when helical turbulence is induced at the scale fe,, 
large scale quantities peaked at a smaller ki (see also Maron 
& Blackman 2002). Therefore we write the different terms of 
equations Clll, 113 and (|3 as: (a.b>,,o/ = /if, <Ao.Bo),.„/ = 
//f, <u.b),,o/ = h^, <b.(Vxb))o = k^hf, (u.(Vxu))o = h", 
<u.V2b>o = -klh^, (b.V^u),) = -klh^, <(V X b) . (V X u))o = 
klh^, <(Vxb).V2b)() = -klhf, <b.(VxV2b))o = -k^hf, 
(V ■ (uxb))o - ±ks\sQ\, V^Bo - -k^fio. Besides we write 
(u.u)o - 2e". We see that in this case the the interaction be- 
tween u and b is described by the small scale cross-helicity 
(u ■ b)o = h'^ ■ For Bq, as it is force-fee, we have |Bo| = 
^V2|;jM|i/2^ and (V X B)) .B,, = /t^/if . When we replace these 
expressions in eqs. and J12I I and J14> we obtain: 



dt 



\[kihf-h^^€„kih^]kfm^ 

2 h'^ 
-ie"±iuksM)kY'^^ 



/2 



[h^-{uk]hf\k][^\h^\'i^-^y^.. 



3 



d 

dt ^ 
and 

dt ' 



h'l^2kf^^^\hf\'l^-2r^kiht 



2,M 



^-2kf^^„m'l'-2r^kih 



(15) 
(16) 

(17) 



h!^ is not an ideal invariant in Hall-MHD, as can be seen from 
its evolution equation. It is obtained by deriving (u ■ b),io/ with 
respect to time, and using eqs. and and reads 



d{vL ■ b) 



vol 



dt 



= -f;/<(V X U) X (V X b))o • Bo 

- ^H(bx(Vxu))o-(VxBo) 
+ (y + 77)((Vxu)-(Vxb))™,. 



(18) 



In order to close our equation system, we could try to make 
in eq. ( I18> the same approximations used in eq. M \\ . \\2\ 
and (I14> . However they would produce expressions for which 
new equations should be deduced. These new equations in turn 
would produce new terms and so on, thus resulting in a system 
difficult to integrate and hard to interpret physically. Therefore 
we shall proceed as follows. The presence of Hall effect im- 
plies that the magnetic field must satisfy V x b 9^ b, i.e. it 
cannot be force-free. Assuming V x b oc b means two things: 
on one side that we are in the standard case (i.e., without Hall 
effect), and on the other, that the small scale magnetic field is 
in an equilibrium state (i.e., no Lorentz force is excerted on 
the stochastic electric currents). Therefore to use V x b oc k^b 
in eq. jl8t means to consider a leading order in a perturba- 
tive expansion of the different terms of eq. M5\ . around the 
Hall-free state, but it does not mean that we are expanding 



around a force-free state, as the approximation is used only 
m eq. ( fTst . There remains the issue of the sign. As we shall 
be interested in a situation in which small scale magnetic he- 
licity grows to negative values, we choose V x b ^ -ks\>, to 
guarantee that condition. For the factor V x v, we shall assume 
maximal negative helicity and thus write it as V x v ^ -ks\. 
We then write the first and second terms in the r.h.s. of eq. 
(US) as ((V X u) X (V X b) ■ Bq),,,,, ^ A:^(u x b)™/ ■ Bq and 
(b X (V X u) ■ (V X Bo)>„o/ - ^5 <b X u)™; ■ (V X Bo). Under this 
approximation, we note that the second term becomes smaller 
than the first one by a factor ki/ks, and that therefore can be 
discarded if ki/ks «; 1. From the remaining expression, we 
see that we would also need the evolution equation for eo ■ Bo. 
However as we are treating Hall effect as a perturbation, we 
make a negligible error if we use sh ■ Bq instead of eo ■ Bo. We 
use the same reasoning to write \sh\ instead of |eol in the second 
term between brackets in the rh.s. of eq. il5\ . We are then left 
with the following equation for the cross helicity: 



dt 



-eHkie"„k]!^\H''\''' 



2lC 



{v + Tf)kiH 



(19) 



and our equation system consists of eqs. M5\ . M6\ . Mil and 
(I19> . In order to numerically integrate it and to correctly de- 
vise the perturbative treatement of the Hall effect, we need to 
make the equations nondimensional. We then define the fol- 
lowing dimensionless quantities: r = ukit, - h'^kiju^, 
g'^ = hfkju^, = h'=lu^, g" = h"l[kLU^), ej!^ = eJ!^/M2, 
Ah = (nkL, ^ = (Ik^^u, f" = e"/u\ Rm = uj {k^rj), r = kslk^^. 
This scheme of normalization is similar to the one of Blackman 



& Field J2002> . except that we use ki instead of ks . Besides 
we shall consider magnetic Prandtl number v/77 = 1 and thus 
ixi + v) - 2rj. When we replace these quantities in eqs. (I15> and 
(I17> -( I19> we obtain the following system: 



dT 



^2|^M|l/2 



- AH[g'-AHr'g^]\GT'-^iQl.. 



do" 



dr " Rur 



^ = -2eLiG«r/2-|^^«, 

dr " Rm 

dT ^" RM^ 



(20) 



(21) 



(22) 



(23) 



In order to have a simple picture of how en evolves, we can 
reason as follows. For high Rm assume a prescribed, negative 
value for g" in the a term of eq. j20> (i.e, the first term between 
square brackets), and that the term A^g'^ is negligible as well as 
the initial value of g'^ . In that situation G*' will initially grow 
toward positive values (due to the first term of eq. 1211 ) and g'^ 
toward negative values (due to the first term of eq. 1221 '). This 
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Fig. 1. Electromotive force for ^ - I, i.e. strong three-point 
correlations. 

will cause the a term to go to zero at a certain instant, and hence 
to the end of the kinetic regime, i.e., the period during which 
the growth of Bq is exponential. The presence of the Ang'' term 
drastically modifies this scenario: were this term negative, then 
it would take the g'^ term a shorter time to cancel the other two 
terms, i.e., we would have a shorter kinetic phase. Were it pos- 
itive, then the opposite situation would occur: the kinetic phase 
would last longer'*. This simple picture is even more modified 
by the fact that now the /3 term in eq. J20> (the second term be- 
tween brackets) is not possitive definite, a fact that could act in 
favour or against of the two situations described above. We can 
conclude that the operation of a Hall-MHD dynamo is far more 
subtle and complicated than the standard MHD one. 

4.2. Discussion 

To integrate the system we used a 4''' order Runge-Kutta 
method with variable step and considered the following val- 
ues for the different parameters that enter in the equations: 
r = 5, = -5 (this value is equivalent to the g" = -1 of 
Blackman & Field I20D2I with the normalization they used), 
/« = 1, 7?^ = 2000 and = and 0.4. The second value 
of Afj corresponds to a Hall length almost twice the turbulent 
scale, but shorter than the coherence large scale. The high value 
of Rm is easily found in astrophysical environments. For the 
three-point correlations we considered two cases: ^ = 1 (strong 
correlations) and { = 2/Rm (weak correlations). As initial con- 
ditions we assumed = 0, g^ - 0.001 - Gq (i.e., an 
initial state with small magnetic helicity. Other initial condi- 
tions do not give qualitative different results) and g'- - 0. We 
also considered the two possible signs in the /3 term, namely 
= /« ± AnrlQll (see eq. (|20j). 

In Fig^we plotted Q^l^ as a function of time, for strong non- 
linearities, i.e. ^ = 1. We see that for Ah + 0, g'j^ is damped 
faster than for Au - 0. For this process is in turn slightly 
stronger than for The saturation value however is the same 
for all cases, i.e. is not modified by Hall effect. The rise of 
the first oscilation in the transitory regime corresponds to the 
kinetic regime, in which the large scale magnetic field would 

From the negative sign of the first term in eq. <23> we see that in 
this case, the first situation will occur, i.e. an earlier quenching of the 
dynamo 



Fig. 2. Magnetic helicities for ^ = 1, i.e. strong three-point 
correlations. Upper curves correspond to large scale MH, and 
lower curves to small scale MH. 

grow exponentially. We see that the instant at which this rise 
stops is slightly smaller than the one at which stops the stan- 
dard MHD curve, while the amplitudes of the curves are sub- 
stantially smaller This instant is independet of either yS+ or 
This behaviour can be interpreted as that the Hall dynamo is 
less efficient than its standard MHD counterpart to generate 
large scale fields. In Fig.|2we plotted small scale magnetic he- 
licity g'^ (lower curves) and large scale magnetic helicity G'^ 
(upper curves), also for ^ = 1. We see that for An + ^ the satu- 
ration value of G'^ is substantially smaller than for An - 0, and 
is the same for both possible yS's. This means that the inverse 
cascade of magnetic helicity is suppressed compared to stan- 
dard MHD, for the considered parameters. Consistently with 
Fig.n we see again that the rise of the first peak takes slightly 
less time for An +^ than for An - 0, and the amplitudes in the 
former case are much smaller than in the latter case. 

In Fig. 13 we plotted the e.m.f. qJ^q for f - 2IRm, i.e. 
weak non-linearities. The smaller amplitudes of the Hall-MHD 
curves means that the e.m.f. is more quenched than for standard 
MHD, as in the case of strong non-linearities. Again in this case 
the quenching due to /?+ is slightly stronger than the one pro- 
duced by In this case the action of Hall effect in the e.m.f. 
is manifested for all times, as no saturation value is attained. 
Again here the rise of the first peak coiTesponds to the kinetic 
regime, and similar features as for ^ = 1 are found: durantion 
slightly shorter and amplitude significantly smaller, showing 
that in this case again the Hall-dynamo is less efficient than the 
standard MHD one. In Fig.0]we plotted the small scale mag- 
netic helicity g'^ (lower curves) and the large scale magnetic 
helicity G*' (upper curves). We can appreciate more clearly the 
quenching produced by Hall effect, as the amplitude of the Hall 
dynamo is about 5 times shorter than the standard MHD. The 
comments about the durantion of the kinetic phase are the same 
as for the ^ = 1 case. 

In all figures, the main features (oscillations for the ^ - 
21 Rm and saturation for ( - 1) are determined by the a term, 
i.e. by the interplay between kinetic and current small scale he- 
licites, and for the Hall dynamo also by the coupling between 
small scale magnetic and velocity fields: Given a prescribed 
negative value of g" in the a term of eq. (I20> . then g'^ grows 
negative, thus leading to a cancellation of them and to a con- 
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Fig. 3. Electromotive force for ^ = 2/RM, i.e. weak three-point 
correlations. 

sequent supression of the growth of gjj^g (this is the end of the 
kinetic regime). In the Hall-dynamo, in the approximation we 
work with, the a term gets an extra term proportional to g*^, 
which in principle can be positive or negative. For the parame- 
ters we considered in this work it is negative and consequently 
reinforces the action of g*^, thus leading to suppression of the 
a term faster than in the non-Hall case. Were it positive, then 
the opposite would occur: it would reinforce -g" and thus it 
would take longer for the (negative) small scale magnetic he- 
licity to catch up with the other (positive) terms in the a term. 
As a result the kinetic regime would last longer. 

The results quoted in this paper, namely a stronger quench- 
ing of the Hall-MHD electromotive force for tiurbuience < 
^Haii < {.system, agree with numerical simulations performed by 
Mininni et al (Mininni, Gomez & Mahaian l2003bt . 

5. Conclusions 

In this paper we studied semianalitically how Hall effect mod- 
ifies the quenching process of the electromotive force in Mean 
Field Dynamo theory. We used a dynamical closure scheme 
named minimal r approximation (Blackman & Field 2002 
Brandenburg & Subramanian .2005J . that takes into account 
the back-reaction of the small scale fields generated by the 
turbulence. As we considered helical turbulence, those small 
scale fields are not to be considered as the result of a small 
scale dynamo, but as a waste product, that nevertheless strongly 
quenches the e.m.f. 

Considering force-free large scale magnetic fields, we 
found that Hall effect modifies the evolution equation of the 
e.m.f. in several ways: the main driving term, the so-called a 
term, proportional to |B()| now depends on the electronic veloc- 
ity Ue instead of the fluid velocity. Besides there appears a third 
term, explicitly dependent on the Hall parameter {h, that cou- 
ples the magnetic field to this electronic velocity. The diffusive 
term, also known as /3 term, proportional to |V x BqI also de- 
pends on the electronic velocity but besides it acquired a new 
term that couples the small scale magnetic field to the elec- 
tronic velocity, and that renders it not positive definite. A neg- 
ative value of this coefficient represents non local transfer from 
small scale turbulent fields to the large scale magnetic field. 
Finally there appears a new term proportional to V^Bp. In our 




T 



Fig. 4. Magnetic helicities for ^ = 2/RM, i.e. weak three-point 
correlations. 



case this term plays no significant role because it is substan- 
tially smaller than the others due to the perturbative scheme we 
use. 

To give a concrete numerical example, we considered Hall 
effect as a perturbation of characteristic scale larger than the 
turbulent scale, but shorter than the large scale magnetic field. 
This situation can be found in e.g. accretion disks and in the 
early universe plasma (Sano & Stone l2002l Tajima et al. ll992> . 
After implementing the two scale approximation, we numeri- 
cally integrated the resulting evolution equations for the e.m.f., 
the large and small scale magnetic helicities, and the cross- 
helicity which is the quantity that in this approximation mim- 
ics the coupling between the small scale velocity and magnetic 
fields. The overall effect is that in the presence of Hall effect, 
the e.m.f. is more strongly quenched than in the case of stan- 
dard MHD dynamo. This fact is acompanied by a damping in 
the inverse cascade of magnetic helicity. 

As the main scope of this paper is to understand concep- 
tually how Hall effect acts on a MFD, we did not apply our 
results to a concrete astrophysical object. We leave this issue 
for future work, after a deeper understanding of the mecha- 
nism is attained. Besides this point, this work can be improved 
in several aspects. A very important one is to extend the per- 
turbative expansion to higher orders (or even to device a non- 
perturbative scheme) in order to attain other regimes, where dy- 
namo action might be enhanced by Hall effect. This extension 
might also mean that we should abandon the hypothesis that Bq 
is force-free. Next, to work with more general boundary con- 
ditions, that permit to relax the conservation of total magnetic 
helicity. Finally, a more detailed study of the three-point corre- 
lations is in order, as well as to consider non-helical turbulence. 
We are working on some of these topics at present. 
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Appendix A: Evolution Equation for tlie 
Electromotive Force 

The evolution equation for is obtained by calculating 



— = (-xb)o-.<ux-)o 

(^lL^xb)o.((Vxb)xf)o , (A.1) 
at at 



and the corresponding equation for e^^ - fiH'(Bo/|Bo|) by doing 
5^/5? = (5£H/5f)(Bo/|Bo|) + fiH[5(Bo/|Bo|)/5f]. Replacing 
eqs. (jSJl and considering the development of operator P to 
first order in the terms linear in Bo, as is done in Gruzinov & 
Diamond ( 1995 1, and Blackman & Field J2002t . and assuming 
homogeneous and isotropic turbulence, we obtain: 

^ = i {<(V X b) ■ b)o - <U ■ (V X u))o - ^H<U ■ V2b)o 
+ ^H<b • V2u)o + ^H<(V X b) ■ (V X u))o 

- 4(b ■ [V X v2b])o + 4((V X b) ■ v2b)o} Bo 

- \ [{v\ + (ni^ ■ (u X b))o] (V X Bo) 

+ i [(u ■ b)o - ^H<(V X b) ■ b)o] V^Bo 

+ 77 [(u, X v2b)o - ^H<(V X V^b) X b)o] 

+ v(V2uxb)o+T, (A.2) 

where by T we denote the non linear terms, i.e.; 

T = (u X [V X (u X b)])o + <P([u X (V X u)] x b))o 

+ <P([(V X b) X b] X b))o - £h{u X {V X [(V X b) X b]!)o 

- {niC^ X b) X [V X (u X b)])o - ^H<{V X [V X (u X b)]) X b)o 



+ 4<(Vxb)x{Vx[(Vxb)xb]!)o 
+ 4<[Vx{Vx[(Vxb)xb])]xb)o. 



(A.3) 



Recalling that = u - ^//V x b, we can write eq. ( IA.2> in a 
form that shows explicitly that now it is the electronic flow the 
driver of the dynamo: 



^ = 1 {<(V X b) ■ b)o - <u, ■ (V X u,))o + {Hih ■ v2u,)o} Bo 



3 

i [(u ■ u,)o + {nib ■ (V X u))o] (V x Bo) 



i<u, ■ b)oV2Bo + 7, [(u, X V2b)o - ^H<(V X V^b) X b)ol 

(A.4) 



3 

+ v(V^uxb)o + T 



